Astronomy & Astrophysics manuscript no. born'rev'ult 


February 2, 2008 


(DOI: will be inserted by hand later) 





The formation and evolution of hydrogen-deficient post- AG B 

white dwarfs: the emerging chemical profile and the 
expectations for the PG1159-DB-DQ evolutionary connection 

L. G. Althaus 1 *, A. M. Serenelli 2 , J. A. Panei 3 **, A. H. Corsico 3 ***, E. Garcia-Berro 1 ' 5 , 

and C. G. Scoccola 3 ^ 

in 
o 

1 Departament de Fi'sica Aplicada, Universitat Politecnica de Catalunya, Av. del Canal Olhnpic, s/n, 08860, 
£SJ ' Castelldefels, Barcelona, Spain 

2 Institute for Advanced Study, School of Natural Sciences, Einstein Drive, Princeton, NJ, 08540, USA 

3 Facultad de Ciencias Astronomicas y Geofi'sicas, Universidad Nacional de La Plata, Paseo del Bosque S/N, 
(B1900FWA) La Plata, Argentina. 

4 Instituto de Astroh'sica La Plata, IALP, CONICET-UNLP 
' 5 Institut d'Estudis Espacials de Catalunya, Ed. Nexus, c/Gran Capita 2, 08034, Barcelona, Spain. 

e-mail: althaus , panei, acorsico, cscoccolaOf caglp.unlp.edu. ar, aldos@ias.edu, garcia@fa.upc.es 

Received; accepted 

in ' 



o 
o 

(N 



Abstract. This paper is designed to explore the formation and evolution of hydrogen-deficient post-AGB white 
dwarfs. To this end, we compute the complete evolution of an initially 2.7 Mq star from the zero-age main 
sequence through the thermally pulsing and mass-loss phases to the white dwarf stage. Particular attention is given 
ly-^ ' to the chemical abundance changes during the whole evolution. A time-dependent scheme for the simultaneous 

treatment of abundance changes caused by nuclear reactions, diffusive overshooting, salt fingers and convection 
is considered. We employed the double-diffusive mixing-length theory of convection for fluids with composition 
gradients (Grossman & Taam 1996). The study can therefore be considered as a test of its performance in low-mass 
stars. Also, time-dependent element diffusion for multicomponent gases is taken into account during the white 
dwarf evolution. The evolutionary stages corresponding to the last helium thermal pulse on the early white-dwarf 
cooling branch and the following born-again episode are carefully explored. Relevant aspects for PG1159 stars and 
DB white dwarf evolution are studied in the frame of these new evolutionary models that take into account the 
history of the white dwarf progenitor. The scope of the calculations is extended to the domain of the helium-rich, 
carbon-contaminated DQ white dwarfs with the aim of exploring the plausibility of the evolutionary connection 
PG1159-DB-DQ. In this regard, the implications for the double-layered chemical structure in pulsating DB white 
r> ' dwarfs is investigated. Another aspect of the investigation concerned the consequences of mass-loss episodes during 

the PG1159 stage for the chemical stratification of the outer layer of DB and DQ white dwarfs. 

Key words, stars: evolution — stars: abundances — stars: AGB stars: interiors — stars: white dwarfs — stars: 
oscillations 
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1. Introduction 

White-dwarf stars with helium-rich atmospheres, com- 
monly referred to as DB white dwarfs, comprise about 
20% of the total white dwarf population. Most of these 
stars are widely believed to be the result of a born-again 
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episode, which is considered to be the most promising 
scenario to explain the existence of hydrogen-deficient 
post-asymptotic-giant-branch (post-AGB) stars — see, 
for instance, Fujimoto (1977), Schonberner (1979) and 
Iben et al. (1983) for earlier references. Other possible 
channels that could lead to the formation of DB white 
dwarfs would involve the class of helium-rich, supergiant 
R Coronae Borealis (RCrB) stars, probably linked to the 
hot hydrogen-deficient O(He) stars, and the hydrogen- 
poor, carbon-normal star evolution directly from the ex- 
tended horizontal branch (AGB manque stars such as 
SdB stars) into the white dwarf state (Schonberner 1996; 
Werner 2001). Within the born-again scenario, a very late 
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helium-shell flash is experienced by a white-dwarf rem- 
nant during its early cooling phase after hydrogen burn- 
ing has almost ceased. At the beginning of this thermal 
pulse, most of the residual hydrogen envelope is engulfed 
by the helium-flash convection zone and completely burnt. 
The star is then forced to evolve rapidly back to the AGB 
and finally into the central star of a planetary nebula at 
high effective temperatures (T e g) but now as a hydrogen- 
deficient, quiescent helium-burning object. Such objects 
are expected to exhibit surface layers that are substan- 
tially enriched with the products of helium burning, par- 
ticularly carbon. 

Important observed examples of these hydrogen-defi- 
cient post-AGB stars are the very hot PG1159 and their 
probable progenitors, the Wolf-Rayet type central stars of 
planetary nebulae having spectral type [WC] 1 (Koesterke 
& Hamann 1997; Dreizler & Heber 1998; Werner 2001). 
Indeed, spectroscopic analyses have revealed that most 
of these post-AGB stars are characterized by hydrogen- 
deficient and helium-, carbon- and oxygen-rich surface 
abundances. In particular, the appreciable abundance of 
oxygen in the atmospheres of these stars has been suc- 
cessfully explained by Herwig et al. (1999) on the basis of 
evolutionary calculations of the born-again scenario that 
incorporate convective overshoot. 

Strong observational evidence suggests that PG1159 
stars are the direct predecessors of the majority of helium- 
rich DO stars which are the hot and immediate progenitors 
of DB white dwarfs (Dreizler & Werner 1996; Dreizler & 
Heber 1998). Theoretical evidence for the existence of an 
evolutionary link between PG1159 stars and most of the 
DO white dwarfs has been presented by Unglaub & Bues 
(2000) on the basis of diffusion calculations with mass loss 
for hot white dwarfs. In addition, evolutionary calcula- 
tions taking into account time-dependent element diffu- 
sion (Dchncr & Kawaler 1995; Gautschy & Althaus 2002) 
have shown that, as a result of gravitational settling of 
carbon and oxygen, the PG1159-likc initial chemical strat- 
ification of a pre-white dwarf evolves into a superficially 
helium dominated double-layered chemical structure when 
the domain of the pulsating DBs is reached. In fact, two 
different chemical transition zones would characterize the 
envelope of the PG1159 descendants: a still uniform inter- 
shell region which is rich in helium, carbon and oxygen; 
the relics of the short-lived mixing episode that occurred 
during the last helium thermal pulse, and an overlying 
pure helium mantle which thickens as cooling proceeds. 
These works clearly foster the plausibility of an evolution- 
ary connection between most of the PG1159 and DB stars. 

The shape of the outer layer chemical profile is a mat- 
ter of the utmost importance as far as the pulsational 
properties of DB white dwarfs is concerned. In fact, the 
presence of a diffusion-induced double-layered chemical 
structure has been shown by Fontaine & Brassard (2002) 
to have strong implications for the theoretical pulsational 
spectrum of these stars. According to these authors, as- 

1 Not to be confused with massive Wolf-Rayet stars. 



teroseismological inference about core composition of DBs 
and the 12 C(a, 7) 16 reaction rate based on single-layered 
DB models (with a pure helium envelope atop a carbon- 
oxygen core) should be taken with a pinch of salt. More re- 
cently, DB asteroseismological fittings incorporating both 
the double-layered envelope feature expected from time- 
dependent diffusion calculations and adjustable carbon- 
oxygen cores have been presented by Metcalfe et al. (2003) 
for a wide range of helium contents and stellar masses. 
Despite their models yielding significantly better fits to 
the observations, the derived stellar parameters for some 
fittings lead them to conclude that double-layered models 
with adjustable carbon-oxygen cores may not be entirely 
appropriate to explain the observations. 

The possibility that PG1159 stars could eventually 
evolve into DB white dwarfs that are characterized by en- 
velopes with a single composition transition zone has re- 
cently been explored by Althaus & Corsico (2004) . Indeed, 
on the basis of evolutionary calculations that incorpo- 
rate time-dependent element diffusion, Althaus & Corsico 
(2004) have shown that if the helium content in PG1159 
stars is smaller than w 10~ 3 M_, the double-layered struc- 
ture is expected to become single-layered by the time 
evolution has proceeded to the domain of the variable 
DBs. Although the small quoted value for the helium 
content is difficult to reconcile with evolutionary calcu- 
lations for the formation of hydrogen-deficient post-AGB, 
which predict the total helium mass left in the star to 
be about M Ho « 10" 2 M» (Herwig et al. 1999), the ex- 
istence of PG1159 stars with low helium content cannot 
be discarded. In fact, mass-loss rates ranging from 10 
to 10 _8 MQ/yr are observed in many luminous PG1159 
stars. In addition, tentative evidence for the persistence 
of mass-loss rates within the range 10~ 7 to 10 -10 MQ/yr 
down to the domain of hot helium-rich white dwarfs has 
been presented (Werner 2001). The existence of such mass- 
loss rates would imply that most of the helium-rich en- 
velope of DB progenitors could be substantially reduced 
during the time interval mass loss would be operative. It 
is worth mentioning that the existence of PG1159 stars 
with a helium content as low as 10 -3 Mq has been sug- 
gested by asteroseismology in at least one of these stars 
with a stellar mass of 0.6 Mq (Kawaler & Bradley 1994), 
thus implying the occurrence of modest mass-loss during 
the evolution to the PG1159 phase. 

In this work we explore some relevant aspects for the 
evolution of DB white dwarfs on the basis of new evolu- 
tionary calculations that account for a complete and self- 
consistent treatment of the evolutionary stages prior to 
the white dwarf formation. We concentrate on DB white 
dwarfs resulting from a born-again episode. Specifically, 
we follow the evolution of an initially 2.7 Mq star from 
the zero-age main sequence through the thermally puls- 
ing and mass-loss phases on the AGB to the white dwarf 
regime. The evolutionary stages corresponding to the al- 
most complete burning of protons following the occur- 
rence of the very late thermal pulse and the ensuing born- 
again episode are carefully explored. Attention is paid to 
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the abundance changes along the whole evolution, which 
are described by means of a time-dependent scheme for 
the simultaneous treatment of nuclear evolution and mix- 
ing processes due to convection, salt fingers and diffusive 
overshoot. We emphasize in particular the role of time- 
dependent element diffusion in the chemical abundance 
distribution in the white-dwarf regime. We also investi- 
gate the influence of mass-loss episodes during the PG1159 
and DO phases for the chemical stratification of pulsating 
DBs. Finally, we extend the scope of our calculations to 
the domain of the helium-rich carbon-contaminated DQ 
white dwarfs, the supposed cooler descendants of DBs. 
The plausibility of the evolutionary connection PG1159- 
DB-DQ (Fontaine & Brassard 2002) is assessed in the 
framework of our new evolutionary models. 

As far as we are aware, this is the first time that the 
evolution of hydrogen-deficient white dwarfs is performed 
consistently on the basis of a complete and detailed treat- 
ment of the physical processes that lead to the formation 
of such stars. At this point it is important to note that 
we think that a re-examination of the pulsational proper- 
ties of variable DB white dwarfs deserves to be performed 
in the frame of the new evolutionary models presented in 
this work. However, being this issue important it is also 
true that this would carry us too far afield. The paper is 
organized as follows. The following section contains de- 
tails on the main physical inputs to the models, particu- 
larly regarding the treatment of the chemical abundance 
changes. In Sect. 3 we present the evolutionary results. 
There, we elaborate on the main aspects of the pre-white 
dwarf evolution, particularly during the born-again phase 
and the attendant chemical changes. In the same section, 
we also describe the results corresponding to the PG1159 
and white-dwarf regimes. In Sect. 4, we discuss the impli- 
cations of our results for the white dwarf evolution as well 
as the role played by mass-loss episodes. Finally, Sect. 5 
is devoted to discuss some concluding remarks. 

2. Input physics and evolutionary sequences 

The calculations presented in this work have been done 
using a Hcnyey type stellar evolution code. In particu- 
lar, we employed the LPCODE evolutionary code that is 
specifically designed to compute the formation and evo- 
lution of white dwarf stars, following a star from the 
main sequence through the thermal pulses and post-AGB 
phases. Except for minor modifications, the code is essen- 
tially that described at length in Althaus et al. (2003) and 
references therein. LPCODE uses OPAL radiative opaci- 
ties (including carbon- and oxygen-rich compositions) for 
different metallicities (Iglesias & Rogers 1996), comple- 
mented, at low temperatures, with the molecular opaci- 
ties from Alexander & Ferguson (1994). Opacities for dif- 
ferent metallicities are required in particular during the 
white dwarf regime, where metallicity gradients induced 
by gravitational settling develop in the envelope of such 
stars. High-density conductive opacities are those of Itoh 
et al. (1994) and the references cited there, whereas neu- 



trino emission rates are those of Itoh (1997) and refer- 
ences therein. The equation of state for the low-density 
regime includes partial ionization for hydrogen and he- 
lium compositions, radiation pressure and ionic contribu- 
tions. For the high-density regime, partially degenerate 
electrons and Coulomb interactions are also considered. 
During the born-again episode, the remnant star develops 
surface layers rich in helium, carbon and oxygen. In that 
case, we employ an ideal equation of state that includes 
partial ionization for any mixture of the three chemical 
species. Finally, for the white dwarf regime, we consid- 
ered an updated version of the equation of state of Magni 
& Mazzitelli (1979). 

The nuclear network employed in LPCODE accounts 
explicitly for the following 16 elements: 1 H, 2 H, 3 He, 4 Hc, 
7 Li, 7 Be, 12 C, 13 C, 14 N, 15 N, 16 0, 17 0, 18 0, 19 F, 20 Ne 
and 22 Ne. In addition, we consider 34 thermonuclear 
reaction rates to describe the hydrogen (proton-proton 
chain and CNO bi-cycle) and helium burning, and carbon 
ignition. Specifically, for hydrogen burning we consider: 
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This set of nuclear reactions allows us to follow in 
detail the main nucleosynthesis occurring during the 
thermally pulsing AGB and born-again phases. Finally, 
the reactions for carbon burning are: 

12 C + 12 C ^ 20 Nc + a 

12 C + 12 C^ 24 Mg + 7 

Nuclear reaction rates are taken from Caughlan & 
Fowler (1988), except for the reactions 15 N(p, 7) 16 0, 
15 N(p,a) 12 C, 18 0(p,a) 15 N, 18 0(p, 7 ) 19 F, 12 C(a, 7 ) 16 0, 
16 O(a, 7 ) 20 Ne, 13 C(a,n) 16 0, 18 0(a, 7 ) 22 Ne, 

22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg, which are taken 
from Angulo et al. (1999). In particular, it is important to 
realize for the purposes of our work that the 12 C(a, 7) 16 
reaction rate given by Angulo et al. (1999) is about twice 
as large as that of Caughlan & Fowler (1988). 

An important aspect of the present study is the treat- 
ment of the abundance changes throughout all the dif- 
ferent evolutionary phases. In particular, during some 
short-lived phases of the evolution, for instance during 
the born-again episode and also the thermally pulsing 
AGB phase, the nuclear time-scale of some reactions be- 
comes comparable to the convective mixing time-scale. In 
this case the instantaneous mixing approximation turns 
out to be completely inadequate for addressing chemical 
mixing in convective regions. A more physically sound 
chemical evolution scheme than instantaneous mixing is 
therefore required. In this work, we have considered a 
time-dependent scheme for the simultaneous treatment of 
chemical changes caused by nuclear burning and mixing 
processes. Specifically, abundance changes are described 
by the set of equations 

fdY\ fdY\ d r. 2 s2 dY 1 

with Y being the vector containing the number fraction of 
all considered elements. Details about the numerical pro- 
cedure for solving this equation can be found in Althaus et 
al. (2003). Here, mixing due to convection, salt fingers and 
overshoot is treated as a diffusion process (second term of 
Eq. ^) . The efficiency of convective and salt-finger mixing 
is described by appropriate diffusion coefficients D which 
are specified by our treatment of convection. In particular, 
we considered the extended mixing-length theory (MLT) 
of convection for fluids with composition gradients devel- 
oped by Grossman et al. (1993) in its local approximation 
as given by Grossman & Taam (1996). These authors have 
developed the non-linear MLT of double-diffusive convec- 
tion 2 that applies in convective, semiconvective and salt- 
finger instability regimes. In this work, the mixing length 
free parameter a is taken to be 1.5. To our knowledge, our 
work constitutes the first application of the extended MLT 
to full evolutionary calculation of low-mass stars. The first 

2 The effects of both thermal and composition gradients de- 
termine the stability of the fluid. 



term of Eq. gives the abundance changes due to ther- 
monuclear reactions. This term is linearized following the 
implicit scheme developed by Arnett & Truran (1969). In 
our treatment of abundance changes, nuclear evolution is 
fully coupled to the current composition change due to the 
various mixing processes — see Althaus et al. (2003) for 
additional details. In LPCODE, abundance changes are 
performed for the 16 elements described above after the 
convergence of each stellar model (and not during itera- 
tions). 

Mixing episodes beyond what is predicted by the 
Schwarzschild criterion for convective stability strongly 
affect the inner chemical profile of white dwarf progeni- 
tors. The occurrence of such mixing episodes, particularly 
core overshooting and/or semiconvection, is suggested by 
both theoretical and observational evidence. Recently, 
Straniero et al. (2003) have presented an assessment of the 
inner chemical abundances in a 3 Mq model star resulting 
from different mixing processes occurring during the late 
stage of core helium burning phase. In particular, they 
conclude that models which incorporate semiconvection 
or a moderate overshoot applied to core and convective 
shells, predict a sharp variation of the chemical composi- 
tion in the carbon-oxygen core. On the other hand, the 
presence of overshooting below the convective envelope 
during the thermal pulses has been shown by Herwig et 
al. (1997) to yield third dredge- up and carbon-rich AGB 
stars for relatively low initial mass progenitors — see also 
Ventura et al. (1999) and Mazzitelli et al. (1999). In addi- 
tion, overshooting below the helium-flash convection zone 
during the thermally pulsing AGB phase gives rise to in- 
tershell abundances in agreement with abundance deter- 
minations in hydrogen-deficient post-AGB remnants such 
as PG1159 stars (Herwig et al. 1999; Herwig 2000). In 
view of these considerations, we have allowed for some 
mechanical overshoot in our work. In particular, we have 
included time-dependent overshoot mixing during all evo- 
lutionary stages. Our scheme for the changes in the abun- 
dances allows for a self-consistent treatment of diffusive 
overshooting in the presence of nuclear burning. We have 
considered exponentially decaying diffusive overshooting 
above and below any formally convective region, includ- 
ing the convective core (main sequence and central helium 
burning phases), the external convective envelope and the 
short-lived helium-flash convection zone which develops 
during the thermal pulses and born-again episode. In par- 
ticular, the expression for the diffusion coefficient in over- 
shoot regions is Z?ov = Dq cxp(— 2z/H v ) where Dq is the 
diffusion coefficient at the edge of the convection zone, z 
is the radial distance from the boundary of the convection 
zone, H v = fHp, where the free parameter / is a mea- 
sure of the extent of the overshoot region, and Hp is the 
pressure scale height at the convective boundary. In this 
study we have adopted / = 0.015, which accounts for the 
observed width of the main sequence and abundances in 
hydrogen-deficient, post-AGB objects (Herwig et al. 1997, 
1999; Herwig 2000; Mazzitelli et al. 1999). We would like 
to mention here that the breathing pulse instability oc- 
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curring towards the end of core helium burning has been 
suppressed — see Straniero et al. (2003) for a recent dis- 
cussion of this point. 

Finally, the evolution of the chemical abundance dis- 
tribution caused by diffusion processes during the white 
dwarf regime has also been taken into account in this 
work. Our time-dependent element diffusion treatment, 
based on the formulation for multicomponent gases pre- 
sented by Burgers (1969), considers gravitational settling, 
chemical and thermal diffusion but not radiative levita- 
tion, which has been neglected, for the nuclear species 
1 H, 3 He, 4 He, 12 C, 13 C, 14 N, 16 and 22 Ne (Althaus et al. 
2003). Diffusion velocities are evaluated at each evolution- 
ary step. We stress that during the white dwarf regime, 
the metal mass fraction Z in the envelope is not assumed 
to be fixed, instead it is specified consistently according 
to the prediction of element diffusion. 

Given the considerable load of computing time de- 
manded by our self-consistent solution of nuclear evolution 
and time-dependent mixing, we have limited ourselves to 
examining only one case for the evolution for the pro- 
genitor star. Specifically, we compute the evolution of an 
initially 2.7 Mq stellar model from the zero-age main se- 
quence all the way from the stages of hydrogen and he- 
lium burning in the core up to the tip of the AGB where 
helium thermal pulses occur. A solar-like initial compo- 
sition (Y,Z) = (0.275,0.02) has been adopted (Anders 
& Grevesse f989). We have considered mass- loss episodes 
taking place during the stages of core helium burning and 
red giant branch following the usual Reimers formula- 
tion with r?R=l. During the thermally pulsing phase we 
adopt the mass-loss rates from Blocker (f995). After ex- 
periencing 10 thermal pulses, the progenitor departs from 
the AGB and evolves towards high effective temperatures. 
Departure from the AGB has been forced to occur at such 
an advanced phase of the helium shell flash cycle that the 
post-AGB remnant undergoes a final thermal pulse during 
the early white dwarf cooling phase — a very late thermal 
pulse, see Blocker (2001) for a review — where most of the 
residual hydrogen envelope is burnt. As mentioned in the 
introduction, there exists observational evidence that sug- 
gests that some hydrogen-deficient post-AGB stars could 
experience strong stellar winds, which could reduce the 
helium content in the star considerably. To assess the in- 
fluence of such mass- loss episodes for the further evolution 
we have considered extreme mass-loss rates of 10~ 7 and 
10~ 8 M Q /yr during the PG1 159 stage. 

The evolutionary calculations have been followed down 
to the domain of the helium-rich carbon-contaminated 
DQ white dwarfs with the aim of exploring the evolu- 
tionary connection between DQ white dwarfs and the 
PGf f 59 stars. The evolutionary sequence from the ZAMS 
to the white dwarf stage comprises rs 70000 stellar models. 
Stellar models are divided into about 1500 mesh points. 
The final mass of the white dwarf remnant is 0.5885 Mq. 
Our models would also be particularly appropriate for pul- 
sational studies of PG1159 and variable DB white dwarfs. 
We report below the main results of our calculations. 



3. Evolutionary results 

3.1. White dwarf progenitor: evolution from the ZAMS 
to the thermally pulsing phase 

In Fig.^we show the complete Hertzsprung- Russell (HR) 
diagram. Our numerical simulation covers all the evolu- 
tionary phases of an initially 2.7 Mq star from the ZAMS 
to the domain of the DQ white dwarfs, including the stages 
corresponding to the helium thermal pulses on the AGB 
and the born-again episode (shown as a dotted line). The 
age (in units of 10 3 yr) counted from the occurrence of the 
last thermal pulse peak on the cooling track is shown at se- 
lected points along the evolutionary track. The total time 
spent in central hydrogen and helium burning amounts to 
6.57 x 10 s yr. After helium is exhausted in the core and 
during the following 1.16 x 10 7 yr, the star evolves towards 
the thermally pulsing phase on the AGB. There, helium 
shell burning becomes unstable and the star experiences 
the well known recurrent thermal instability commonly re- 
ferred to as helium thermal pulses (Schwarzschild & Harm 
1965). After 10 thermal pulses and as a result of strong 
mass loss, the remnant star leaves the AGB and evolves 
towards high T c gs. This takes place when the luminos- 
ity of the star is supported mostly by stationary hydro- 
gen burning. Mass loss decreases the stellar mass from 
2.7 to 0.5885 Mq. Note that the post-AGB remnant un- 
dergoes a further (last) thermal pulse on its early white 
dwarf cooling track shortly after hydrogen burning has 
virtually ceased (Blocker 2001). During this born-again 
episode, most of the residual hydrogen envelope is en- 
gulfed by the deep helium-flash convection zone and is 
completely burnt. Evolution proceeds through these stages 
very fast, since it takes only about 30 yr for the remnant 
to expand from a white dwarf configuration to giant di- 
mensions. Note also that after the born-again episode, the 
now hydrogen-deficient post-AGB remnant experiences a 
second excursion towards lower temperatures (a double- 
loop path) before reaching the domain of the PG1159 stars 
and eventually its terminal white dwarf cooling track. It 
is worth noting that during these stages, the remnant is 
a quiescent helium-burning object that reaches the point 
of maximum T e g (the knee in the HR diagram) in about 
15000 yr after the occurrence of the last thermal pulse. 

The evolution during the thermally pulsing phase is 
documented in Fig. |3 This figure shows the time depen- 
dence of the surface luminosity (L SUT ) and the hydrogen- 
and helium-burning luminosities (Lcno and Lhc respec- 
tively) of our initially 2.7 Mq model star, where the time- 
scale is given in Myr from the main sequence (ZAMS). 
A total of 10 thermal pulses with an interpulse period of 
roughly 1.2 x 10 5 yr have been computed before the rem- 
nant leaves the AGB as a result of the enhanced mass- 
loss. Note that the helium burning rate rises very steeply 
at the peak of each pulse, even during the first pulses. 
In our simulation, departure from the AGB takes place 
at such an advanced stage in the helium shell flash cycle 
that the post-AGB remnant will experience a last helium 
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Fig. 1. Hcrtzsprung-Russell diagram for the complete evolution of our 2.7 Mq stellar model from the ZAMS to the 
white dwarf domain. The star experiences a very last thermal pulse (VLTP) on its early cooling phase after hydrogen 
burning has almost ceased. The following evolutionary stages correspond to the born-again phase and are depicted, 
together with the VLTP stage with dotted lines. The numbers close to the filled dots along the evolutionary track 
correspond to the age (in 10 3 yr) counted from the occurrence of the peak of the last thermal pulse. The domain of the 
PG1159, DO, DB and DQ stars are indicated, as well as the thermal pulse and born-again stages. As a result of mass 
loss episodes, the stellar mass decreases from 2.7 to 0.5885 Mq. After the born-again episode, the hydrogen-deficient 
post-AGB remnant experiences a second excursion towards lower temperatures before reaching eventually its terminal 
white dwarf track. 



thermal pulse at high T e g values — see Fig. Q] and also 
the next section. This situation corresponds to the last 
thermal pulse and it is illustrated separately on the right 
panel of Fig. [21 During the last pulse, hydrogen-burning 
luminosity, due mainly to proton captures by 12 C, reaches 
about 10 8 - 5 Lq. 

The inner carbon, oxygen and helium distribution in 
the progenitor star as a function of the mass coordinate is 
shown in Fig. [3] The upper panel of the figure depicts the 
chemical profile at the onset of the thermally pulsing AGB 
phase, while the bottom panel illustrates the situation for 
the beginning of the last thermal pulse at high T e g. Some 
features of this figure deserve comments. In particular, 



the carbon and oxygen abundance distribution within the 
core is typical of situations in which extra mixing episodes 
beyond the fully convective core during central helium 
burning are allowed. Such extra mixing episodes, partic- 
ularly core overshooting and/or semiconvection, have a 
large influence on the carbon and oxygen distribution in 
the core of white dwarfs (Straniero et al. 2003). In our 
simulation, the central oxygen mass abundance amounts 
to 0.724 when the helium convective core disappears. It is 
worth noting that the shape of the carbon-oxygen profile, 
particularly the sharp variation around M r w 0.33 Mq 
induced by mechanical overshoot, is in agreement with 
that reported in Straniero et al. (2003). The importance 
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Fig. 2. The temporal evolution of surface luminosity and 
hydrogen- and helium-burning luminosities in solar units 
for an initially 2.7 Mq star during the thermally pulsing 
phase on the AGB. After experiencing 10 thermal pulses, 
the progenitor star departs from the AGB with a thermal- 
pulse cycle phase of 0.89 (denoted by a black dot). The last 
thermal pulse shown here takes place already on the white 
dwarf cooling track, where the hydrogen-burning luminos- 
ity exceeds 10 8 Lq. The right panel shows the luminosity 
evolution after the occurrence of the tenth pulse. Age is 
given in 10 6 yr counted from the main sequence. 



of such extra mixing episodes, particularly core overshoot, 
for the pulsational properties of massive ZZ Ceti stars has 
recently been emphasized by Althaus et al. (2003) and 
Corsico et al. (2004) . Note that during the thermally puls- 
ing phase, the mass of the carbon-oxygen core increases 
considerably because the helium burning shell moves out- 
wards. The signatures left by overshoot in the chemical 
profile during this stage is apparent. In particular, the 
overshoot region stretching below the short-lived helium- 
flash convection zone at pulse peak causes the intershell 
region below the almost pure helium buffer to be substan- 
tially enriched in oxygen. Another feature predicted by 
our calculations is the occurrence of the third dredge-up 
during the thermally pulsing phase. We will not discuss 
this aspect here — see Herwig (2000) for details — suffice 
it to say that as a result of envelope overshoot, stretching 
down to the underlying carbon-rich layers, large amounts 
of carbon are dredged up to the surface after each thermal 
pulse. Indeed, the carbon-oxygen ratio increases from « 
0.25 when entering the thermally pulsing phase to « 0.34 
by the time the remnant departs from the AGB. It is im- 
portant to point out at this point that we have checked the 
formation of a carbon star in our simulation by suppress- 
ing mass loss. We find in that case that a carbon-oxygen 
ratio larger than 1 occurs after the occurrence of the 15 th 
pulse. 

Finally, another prediction of our calculations is the 
formation of 13 C- and 14 N- pockets at the base of the he- 
lium buffer after the end of the dredge-up phase. Indeed, 
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Fig. 3. Internal 4 He, 12 C and 16 profiles for the 2.7 M Q 
DB white dwarf progenitor shortly before the occurrence 
of the first thermal pulse on the AGB (upper panel) and 
at the onset of the last thermal pulse at high T e g (bot- 
tom panel). Overshoot episodes leave recognizable features 
both in the inner carbon-oxygen profile and in the inter- 
shell region that are particularly rich in oxygen. See text 
for details. 

during the third dredge-up, diffusive overshoot has formed 
a small region in which hydrogen from the envelope and 
carbon (resulting from helium burning) from the intershell 
region coexist in appreciable abundances. When this re- 
gion heats up enough, hydrogen reignites giving rise to the 
formation of a 14 N-pocket with 14 N abundances by mass 
of about 0.45 (the most abundant element in the pocket) 
and an underlying 13 C-pocket with a maximum 13 C abun- 
dance of 0.06. The mass range over which the 14 N-rich 
region extends amounts only to « 2 x 10~ 6 ,AfQ. During 
the interpulse period, the 13 C-pocket is radiatively burnt 
at relatively low temperatures (about 80 x 10 6 K) before 
the onset of the next pulse via the reaction 13 C(a,n) 16 0, 
the main neutron source reaction in AGB stars. The 
14 N-pocket is engulfed by the helium-flash convection 
zone during the next thermal pulse and burnt via the 
14 N(a, 7 ) 18 F(/3+,j/) 18 0(a,7) 22 Ne chain that converts all 
the abundantly present 14 N to 22 Ne. In agreement with 
Herwig (2000), our calculations show that diffusive over- 
shooting at the base of the convective envelope is a pro- 
cess that naturally leads to the formation of 13 C- and 14 N- 
pockets, which is a fundamental issue regarding the forma- 
tion of heavy elements through the slow neutron capture 
process — see, for instance, Lugaro et al. (2003) for a 
recent discussion. 

3.2. Evolution during the last thermal pulse and born 
again: the burnt of hydrogen envelope 

As a result of the mass-loss episodes, the mass of the hy- 
drogen envelope is reduced to such an extent that the 
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Fig. 4. Hertzsprung-Russell diagram for the evolutionary stages following the onset of the last helium thermal pulse 
that takes place at high effective temperatures (born-again episode) for our post-AGB remnant of 0.5885 Mq. Note 
that the remnant experiences a second excursion (double loop) to the red giant region before evolving into its final 
white dwarf state. The left inset displays the time-dependence of hydrogen- and helium-burning luminosities (solid and 
dotted lines, respectively) expressed in solar units. Selected stages are labeled by letters along the curves. In particular, 
points A and C correspond to the moment of maximum helium- and hydrogen- burning luminosities, respectively. At 
point B, protons begin to be engulfed and burnt by the outward-growing, helium-flash convection zone. Between points 
E and F, hydrogen burning becomes almost extinguished. The hydrogen content remaining in the star amounts to 
1.3 x 10~ 8 Mq. Note the extremely short evolutionary time-scales characterizing this phase. 



thermally pulsing progenitor star leaves the AGB and 
evolves into the planetary nebula regime at large T e g 
values. During this phase of the evolution, helium-shell 
burning increases gradually and when hydrogen burning 
becomes virtually extinct at the beginning of the white 
dwarf cooling branch, the post-AGB remnant experiences 
the last helium-shell flash that gives rise to the short-lived 
born-again episode. The mass of hydrogen that it is left 
in the very outer layers at the start of the last thermal 
pulse amounts to 7 x 10~ 5 Mq. It is worth mentioning 
that very few detailed numerical simulations through this 
complicated regime exist in the literature. In particular, 
evolutionary calculations that include hydrogen and he- 



lium burning combined with time-dependent mixing have 
been performed initially by Iben & MacDonald (1995) and 
by Herwig et al. (1999) for the situation in which diffusive 
overshooting is considered. Also, Lawlor & MacDonald 
(2003) have recently presented a grid of stellar evolution- 
ary calculations for the born-again phenomenon aimed 
to explain the observational characteristics of born-again 
stars. In what follows, we report the main predictions of 
our calculations for the born-again phase, particularly em- 
phasizing on the changes in the abundances of the different 
chemical species that take place during this brief evolu- 
tionary phase. 
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The Hertzsprung-Russell diagram focusing onto the 
last helium thermal pulse and the subsequent born-again 
evolutionary phase is displayed in Fig. 0] Selected evo- 
lutionary stages are labelled by letters along the evolu- 
tionary track. The inset shows the time-dependence of 
hydrogen- (CNO-cycle reactions) and helium-burning lu- 
minosities, again selected evolutionary stages are corre- 
spondingly labeled in this curve. Note that remarkable 
changes in the structure of the star take place on ex- 
tremely short time-scales. For instance, it takes 0.4 yr 
for the star to develop hydrogen-burning luminosities as 
high as 10 s Lq (point C in Fig. 0} after protons begin to 
be engulfed by the outward-growing helium-flash convec- 
tion zone (point B). Most of the hydrogen envelope burn- 
ing occurs between points C and E in about 1 month. 
Specifically, at point E, after 0.9 yr have elapsed from 
the onset of the last helium thermal pulse (point A), the 
mass of the residual hydrogen envelope has reduced to 
10~ 7 Mq. Between points E and F, hydrogen burning be- 
comes gradually extinct and the mass of hydrogen remain- 
ing in the star amounts to 1.3 x 10~ 8 Mq. Approximately 
0.6 yr later, the remnant reaches the point of maximum ef- 
fective temperature at logT c ff = 4.8 for the first time after 
the helium flash. Afterwards, the evolution proceeds into 
the red giant domain somewhat more slowly. In fact, the 
effective temperature decreases to 10000 K over a period 
of about 20 yr, and to 5200K in about 50 yr and the radius 
of the star increases to 30.2 and 125 Rq, respectively. 

It is worth noting that our born-again time-scale of 
20-40 yr is larger than the evolutionary time-scale of the 
born-again Sakurai object (V4334 Sgr), which has evolved 
from the pre-white dwarf stage into a AGB giant star in 
only about 6 yr. The evolutionary calculations of Herwig 
(2001), which are based on the standard MLT, also pre- 
dict too large born-again time-scales (typically 350 yr). 
Herwig has found that the very short born-again evolu- 
tionary time of V4334 Sgr can be reproduced by stellar 
models if the convective mixing efficiency in the helium- 
flash convection zone is reduced by a factor of 100 be- 
low that obtained from the MLT. Although our evolution 
time-scales are larger than observed ones, they are much 
shorter than those obtained by Herwig (2001), but with- 
out invoking an additional reduction in the mixing effi- 
ciency. We remind the reader that in our calculation the 
double-diffusive MLT for fluids with composition gradi- 
ents (Grossman & Taam 1996) has been used 3 . Another 
prediction of our calculations is the occurrence of a double- 
loop in the Hertzsprung-Russell diagram. That is, the star 
reaches red giant dimensions for a second time after the 
onset of the last helium thermal pulse and before finally 
returning to the white dwarf cooling track, a behaviour re- 
ported by Lawlor & MacDonald (2003) and Herwig (2003) 
in the case of low convective mixing efficiency. In partic- 



3 It is clear that the big difference in the evolutionary time 
scales between the both set of calculations cannot be entirely 
attributed to the fact that the mass value of our model is not 
strictly the same as that analyzed by Herwig 2001. 



ular, this second return to the AGB takes about 350 yr 
in our calculations. In the light of these results, we judge 
that a more comprehensive comparison between the stan- 
dard mixing length theory and the double-diffusive mix- 
ing length theory deserves to be done. Such a comparison 
would, however, carry us too far afield and, consequently, 
we postpone it to a forthcoming work. 

During the last helium thermal pulse, profound 
changes in the chemical structure take place as a result of 
the vigorous nuclear processing and mixing episodes. Such 
changes are crucial for the subsequent evolution of the 
star. Because of the very short evolutionary time-scales 
characterizing this phase, a numerical treatment that con- 
sistently couples the equations of nuclear changes with the 
equations for time-dependent mixing processes, like the 
one adopted here, is indeed required for a realistic descrip- 
tion of the abundance changes. A complete coverage of 
the inner chemistry variations that take place during this 
short evolutionary stage is provided in Fig.0 which shows 
the chemical abundance distribution at four selected evo- 
lutionary stages during the last thermal pulse. Specifically, 
the abundances by mass of 1 H, 4 He, 12 C, 13 C, 14 N and 16 
are ploted in terms of the outer mass fraction q. Grey and 
shaded regions mark the domains of convection and over- 
shooting. Panel a shows the chemical stratification at the 
start of helium thermal pulse (point A in figure 4). In the 
outermost layers the chemical composition corresponds to 
that fixed by dredge-up episodes during the AGB phase. 
In the pure helium buffer, the relatively large abundance 
of 14 N reflects the efficiency of hydrogen burning during 
previous evolutionary phases in processing CNO elements 
into 14 N. Because of the large amount of energy resulting 
from helium burning, an outward-growing convective re- 
gion develops. In about 1 yr, the outer edge of this convec- 
tive region reaches the base of the hydrogen-rich envelope. 
As a result, protons begin to be transported downwards 
into hotter and carbon-rich layers, where they are cap- 
tured via the 12 C(p, 7) 13 N reaction. The ensuing vigorous 
hydrogen burning forces the development of an entropy 
barrier. As a result, the original convective region is split 
at logg w —1.6 in two distinct convective shells. The in- 
termediate region is radiative, but also presents a salt- 
finger instability at its top. This situation is illustrated 
by panel b. Hydrogen and helium burning take place at 
the base of such convective layers. Note the large amount 
of 13 C in the convection zone powered by hydrogen burn- 
ing (upper convection zone). In the hotter, helium-flash 
convection zone 13 C (and 14 N) is destroyed by a captures 
and neutrons are released. Note as well that, as previously 
emphasized, the use of a simultaneous treatment of mix- 
ing and burning is indeed required during these stages, as 
reflected by the hydrogen profile. 

Panel c depicts the situation two weeks later. During 
this time, the outer edge of the upper convection zone 
propagates further outwards in mass, mixing downwards 
protons from the original (unprocessed) hydrogen-rich en- 
velope. Note both the persistence of the small salt-finger 
region separating the two convection zones and the in- 
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Fig. 5. Internal abundance distribution of 1 H, 4 He, 12 C, 13 C, 14 N and 16 as a function of the outer mass fraction 
q for the 0.5885 Mq remnant at various selected epochs during the course of the last helium thermal pulse. Panel a 
corresponds to the moment just before the occurrence of the peak helium-burning luminosity (point A in figure 4). The 
situation when the outward-moving, helium-flash convection zone has reached the base of the hydrogen-rich envelope 
(between points B and C in figure 4) is illustrated by panel b. Note that protons are ingested into deeper and hotter 
layers. As a result of the ongoing hydrogen burning a small radiative and salt-finger zone establishes at log q « -1.6, 
which splits the convection zone into two. The situation two weeks later is visualized in panel c. During this time, 
the convection zone powered by hydrogen burning propagates outwards. Finally, panel d illustrates the abundance 
profiles near point F in figure 4. Only traces of hydrogen remain after hydrogen burning. Nuclear energy release due 
to hydrogen and helium burning are also shown in panel b and c. 



crease in the 13 C and 1 N abundances in the upper con- 
vection zone due to hydrogen processing. The surface and 
upper convection zones merge temporarily (point E in Fig. 
0J) , causing the convectively unstable region to extend from 
the hydrogen-burning zone essentially to the very surface 
of the star. Only a very thin unprocessed hydrogen-rich 
envelope remains, which is expected to be diluited by sur- 
face convection when the star returns to the AGB. Finally, 
panel d shows the abundance profiles near point F in Fig. 

Here, hydrogen burning is virtually extinct and the re- 
maining hydrogen mass amounts to only 1.3 x 10~ 8 Mq. 
Note also that traces of hydrogen are present in layers 
even as deep as 3 x 10~ 4 Mq below the stellar surface. 
Interestingly enough, the content of 13 C left in the whole 
star is sizeable and amounts to 4.5 x 1CP 4 Mq, that is 4-5 



orders of magnitude larger than the hydrogen mass of the 
star. 

3.3. Evolution towards the PG1159 state 

After the last thermal pulse, as the star evolves back to 
giant dimensions, either convective dilution or mass loss 
during the quiescent helium-burning phase (Iben et al. 
1983) are expected to erode the tiny layer of original enve- 
lope material, exposing the underlying hydrogen-deficient 
layers. Thus, after the second loop in the HR diagram, the 
remnant evolves to the region of the PG1159 stars with 
a surface chemical composition similar to that shown in 
panel d of Fig. Specifically, in the outer layers, 4 He, 12 C 
and 16 O are by far the dominant species with abundance 
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Fig. 6. Evolution of the internal abundances of 13 C 
and 16 O (dashed and solid lines respectively) across the 
PG1159 domain in terms of the outer mass fraction q for 
the 0.5885 Mq remnant. Thin and thick lines correspond, 
respectively, to a model at T e g = 71000 K before the point 
of maximun effective temperature in the HR diagram and 
T e ff= 125000 K below that point. Note the formation of a 
bump in the oxygen distribution as a result of a captures 
by 13 C. 

by mass ( 4 He, 12 C, 16 0)= (0.306,0.376,0.228). Amongst the 
main remaining constituents are 13 C, 14 N and 22 Ne with 
mass fractions of 4, 1.2 and 2.1 %, respectively. 13 C and 
1 N are present from the outermost layers down to the 
base of the former upper convection zone at 0.01 Mq be- 
low the stellar surface. At deeper layers, these chemical 
species have been already depleted via a captures. The 
large mass fraction of oxygen is an indication of the occur- 
rence of diffusive overshooting during the thermally puls- 
ing AGB phase (see also Herwig et al. 1999). The final 
surface composition of our models is also in line with sur- 
face abundance patterns observed in most hot, hydrogen- 
deficient post AGB stars such as PG1159 stars and central 
stars of planetary nebulae of spectral type WC (Koesterke 
& Hamann 1997; Dreizler & Heber 1998; Werner 2001). 
The remarkable agreement between the 14 N abundance 
predicted by our calculations and the one detected by 
Dreizler & Heber (1998) in five out of nine PG1159 stars 
strongly supports the hypothesis that these stars would 
be AGB descendants that have experienced a born-again 
episode. As we have already seen, mixing and burning of 
protons in the helium-flash convective region are indeed 
required to synthesize 14 N in a carbon- and oxygen-rich 
environment. 

Finally, we note that as the evolution proceeds through 
the domain of the PG1159 stars, the changes in the chem- 
ical composition take place as a result of nuclear burning. 
Fig. is an example of this. Note in particular that in a 
region stretching from 3x 10 -3 to 10~ 2 Mq below the stel- 
lar surface, a bump in the oxygen distribution is built up 
as a result of 13 C burning via the 13 C(a, rt) 16 reaction. 
This is a consequence of the fact that before the hydrogen- 



deficient remnant reaches the point of maximum effective 
temperature in the HR diagram, the temperature at the 
tail of the 13 C distribution exceeds 85 xlO 6 K, which is 
high enough for the 13 C(a, n) 16 reaction to operate. As a 
last remark, we note that during the early PG1159 stage, 
the helium content in the star is reduced from 8.5 x 10 -3 
to 5.2 x 10 -3 Mq as a result of helium burning. The last 
emoted value is the amount of helium with which the star 
enters the white dwarf domain after helium burning be- 
comes extinct. Because we have not invoked any addi- 
tional mass loss after the born-again episode, the amount 
of helium left in the white dwarf (5.2 x 10~ 3 Mq ) should 
be considered as an upper limit for the particular case 
analyzed in this work (see next section for the effect of 
mass-loss episodes during the hot stages of post-AGB evo- 
lution) . 

In Fig. [7J we show as a function of time the luminos- 
ity contributions due to helium burning (Ln c ), neutrino 
losses (-Lneu), surface luminosity (L sur ), and gravothermal 
energy release (L g ) for the hydrogen-deficient 0.5885 Mq 
remnant. The abscissa covers the time span from the pre- 
white dwarf state, before the maximum T c g point in the 
HR diagram, to beyond the domain of the variable DB 
white dwarfs. Luminosities are given in solar units; the age 
is counted from the moment at which the remnant reached 
T e s= 70000 K. Some features are worthy of a comment. 
To begin with, note that at early times it is mostly helium 
burning that contributes to the surface luminosity of the 
star. This is true for the first 15000-20000 yr of evolu- 
tion, shortly after the remnant reaches the maximum T c g 
value, and begins to decrease as the star approaches the 
white dwarf domain. Afterwards, the contribution of he- 
lium burning declines steeply and the evolution of the star 
is dictated essentially by neutrino losses and the release 
of gravothermal energy. At the T e g value characterizing 
PG1159-035 (logT e ff ~ 5.12), gravothermal energy is the 
main energy source of the star. But, for the coolest pul- 
sators in the GW Vir strip neutrino losses exceed photon 
luminosity. This is an important feature since it allows 
to use cool pulsating PG1159 stars to constrain neutrino 
emission processes from measurements of the rate of pe- 
riod change — see O'Brien et al. (1998) and O'Brian & 
Kawaler (2000) for a discussion of this exciting issue. Note 
that neutrino losses constitute the primary cooling mech- 
anism over a period of about 1.5 x 10 7 yr, during which 
the star evolves through the temperature range 90000 K 
;> T e g^ 25000 K. It is worth mentioning at this point that 
the lower temperature limit corresponds to the domain of 
the hot pulsating DB white dwarfs. Thus, these variable 
stars should also be potentially useful to place constraints 
on plasmon neutrino emissivity, as recently emphasized by 
Winget et al. (2004). 

Finally, in Fig.|Hl we show the surface gravity-effective 
temperature diagram for our post born-again 0.5885 Mq 
sequence together with observational data for hydrogen- 
deficient PG1159 and WC stars, as taken from Werner 
et al. (1997) and Werner (2001). In particular, the pul- 
sating nitrogen-rich PG1159 stars analyzed by Dreizler 
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Fig. 7. Time-dependence of the different luminosity contributions (in solar units) for the post born again 0.5885 Mq 
remnant of the evolution of the 2.7 Mq star: surface luminosity, L sur , helium burning luminosity, Lho, neutrino 
losses, L ncu and rate of gravothermal (compressional plus thermal) energy release L g . The evolution of the effective 
temperature (right scale) is also plotted. Time is in years counted from the moment at which the remnant reaches 
T e ff= 70000 K. The vertical lines bracket the domain where neutrino losses exceed photon luminosity. 



& Heber (1998) — namely, PG11 59-035, PG213f+066, 
PG1707+427 and PG0122+200 — are consistent with a 
stellar mass somewhat lower than that characterizing our 
post born-again sequence. 

3.4. White dwarf evolution 

Once helium burning becomes virtually extinct at T e ff« 
100000 K, after 100000 yr of evolution from the last he- 
lium thermal pulse, the hydrogen-deficient remnant set- 
tles upon its terminal cooling track. Here, its chemical 
abundance distribution will be strongly modified by the 
various diffusion processes acting during the white dwarf 
evolution. This can be seen in Fig. |U| which illustrates the 
mass abundances of 4 He, 12 C, 16 0, 13 C, 14 N, and 22 Ne 
for the 0.5885 Mq white dwarf as a function of the outer 
mass fraction at various epochs characterized by values 
of log L/Lq and logT e fj (the corresponding values are 
given in parentheses for each one of the panels). Panel a 
shows the chemical stratification at the start of the cool- 
ing track after the star has reached the point of maximum 
T e ff at high luminosities. In the outermost layers the chem- 
ical composition corresponds essentially to that emerging 



from the mixing and burning events during the last he- 
lium thermal pulse and subsequent born-again episode. 
Rapidly, gravitational settling causes helium to float to 
the surface and heavier elements to sink. In fact, in about 
1.5 x 10 6 yr the star develops a pure helium envelope of 
nearly 6 x 10 -8 Mq (panel b). By the time the domain of 
the variable DBs is reached (after 1.77 x 10 7 yr of white 
dwarf evolution; panel c) gravitationally induced diffusion 
has led to the development of a double-layered chemi- 
cal structure characterized by a pure helium envelope of 
2.3 x 10 -6 Mq atop an intermediate remnant shell rich in 
helium, carbon and oxygen, the relics of the last helium 
thermal pulse (see previous section). About 3.7 x 10 8 yr 
later, the white dwarf reaches the domain of the helium- 
rich, carbon-contaminated DQ white dwarfs. The corre- 
sponding chemical stratification is shown in panel d. Even 
at such an advanced stage, the star is characterized by 
a double-layered structure. In particular, the pure helium 
envelope amounts to 1.9 x 10~ 5 Mq. Note the significant 
carbon enrichment in the surface layers as a result of con- 
vective dredge-up of the carbon diffusive tail by the su- 
perficial helium convection zone (see later in this section) . 
The neutron excess characterizing 13 C (and also 22 Ne) 
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Fig. 8. Evolution of the surface gravity (in cgs units) as 
a function of the effective temperature for the post born- 
again 0.5885 Mq remnant and its comparison with some 
observed hot hydrogen-deficient stars. Filled (open) cir- 
cles denote pulsating (non-pulsating) PG1159 stars, and 
asterisks stand for observed WC stars (Werner et al. 1997; 
Werner 2001). 



partially explains the fact that this element appreciably 
diffuses downwards. It is clear from these figures that diffu- 
sion processes substantially alter the chemical abundance 
distribution in the course of white dwarf evolution. 

The effect of element diffusion on the main chemical 
constituents can best be visualized in Fig. 1101 particularly 
the formation of the double-layered chemical structure. 
The mass abundance of 4 He, 12 C and 16 are shown as 
a function of the outer mass fraction for the 0.5885 Mq 
white dwarf remnant at three evolutionary stages. The 
chemical profiles shortly after the remnant reaches the 
point of maximum T e ff in the HR diagram are represented 
with thin lines having various patterns. Later evolution- 
ary stages around the domain of the DB instability strip 
(log T c ff = 4.41 and 4.28) are represented with normal 
and heavy lines. We stress again that the helium content 
that is left in the star at the start of the cooling branch, 
after helium burning has virtually ceased, amounts to 5.2 
x10 -3 A/q. Because we have not invoked any additional 
mass loss during the hot post-AGB stages or early dur- 
ing the cooling branch, the quoted value for the final he- 
lium mass should be considered as an upper limit for the 
particular case of evolution analyzed here. The diffusion- 
induced double-layered structure at the domain of the pul- 
sating DBs is easily recognizable. Another feature worthy 
of comment is the mixing episode that takes place in the 
region below the intershell zone around logq = — 1. In 



fact, this region is characterized by a inward-decreasing 
mean molecular weight induced by the occurrence of over- 
shooting during the AGB thermally pulsing phase. The 
resulting salt-finger mixing is responsible for the redis- 
tribution of the chemical species in that region, as it is 
apparent from Fig. ITU1 

Finally, we have extended the scope of our evolutionary 
calculations down to the domain of the carbon-enriched 
DQ white dwarfs, thus covering the possible evolution- 
ary connection PG1159-DB-DQ. DQ white dwarfs have 
T eff values below 13000 K (Weidemann & Koester 1995) 
and are characterized by the presence of trace amounts 
of carbon in their atmospheres with abundances by num- 
ber relative to helium, log(nc/riHe)) ranging from —7.3 to 
— 1.5 (MacDonald et al. 1998). The presence of traces of 
carbon observed in these stars is widely believed to re- 
sult from convective dredge-up of the carbon diffusive tail 
by the superficial helium convection zone (Pelletier et al. 
1986; Koester et al. 1982). By the time the DQ domain is 
reached, our white dwarf models have developed a double- 
layered chemical structure with a pure helium mantle of 
1.9 x 10 -5 Mq, which is almost an order of magnitude as 
massive as the helium mantle characterizing our models 
at the beginning of the DB instability strip. However, it is 
less massive than required for heavy elements to be abun- 
dantly dredged-up to the surface by the inwards-growing 
superficial convection zone. This fact is reflected by Fig. 1111 
which shows the surface abundance by number relative to 
4 He of 12 C, 16 O, 13 C and 14 N as a function of the effective 
temperature. As the white dwarf cools down, the base of 
the convection zone moves deeper into the star, with the 
consequent further enrichment with heavy elements of the 
outer layers. Note that 12 C is by far the most abundant 
dredged-up element, with abundances far exceeding the 
low carbon abundances observed in many DQ (see next 
section). Note also that our calculations predict the pres- 
ence of 13 C and 14 N, abundantly created during the last 
helium thermal pulse, as well as 16 in the atmospheres 
of DQ white dwarfs. Finally, 22 Ne has diffused so deep 
into the star (see El panel d) that it is expected not to be 
dredged-up to the surface of these stars. 

For the sake of completeness, we provide in Table 
1 some relevant quantities for our post born-again 
0.5885 Mq sequence. Specifically, we list from left to right 
the effective temperature, the photon luminosity (in so- 
lar units), the age (in years) counted from the moment 
at which the remnant reaches T c s= 10000 K (after the 
occurrence of the double loop at high luminosities), the 
stellar radius (in cm), the surface gravity, and the helium- 
burning and neutrino luminosities (both in solar units). 
The tabulation covers the stages following the end of the 
born-again episode to the domain of the DQ white dwarfs 
at low T e s values. 
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Fig. 9. Abundance by mass of 4 He, 12 C, 16 0, 13 C, 14 N and 22 Ne in terms of the outer mass fraction for the 0.5885 M Q 
remnant at four selected white-dwarf evolutionary stages, characterized by values of log L/Lq and logT e ff (the cor- 
responding values are given in parentheses). Panel a corresponds to the start of the cooling branch, panel c to an 
evolutionary stage during the DB instability strip and panel d to the domain of the carbon-enriched DQ white dwarfs. 
Panel d illustrates the significant carbon enrichment in the surface layers as a result of convective dredge-up of the 
carbon diffusive tail by the superficial helium convection zone. It is clear that element diffusion strongly modifies the 
internal chemical profiles of the hydrogen-deficient white dwarf. 



4. Discussion 

The calculations presented in this work cover the whole 
evolutionary stages involved in the formation of DB white 
dwarfs via the born- again scenario starting on the ZAMS. 
In particular, our calculations have followed the model 
star to very advanced stages of evolution down to the 
domain of the helium-rich carbon-contaminated DQ, the 
supposed cooler descendants of DBs. Hence, the study 
opens the possibility of assessing the evolutionary con- 
nection PG1159-DB-DQ (Fontaine & Brassard 2002) in 
the framework of complete evolutionary calculations that 



take into account a detailed treatment of the physical pro- 
cesses that lead to the formation of hydrogen-deficient 
white dwarfs. In this sense, the results presented here re- 
inforce the conclusions arrived at in Fontaine & Brassard 
(2002) and Althaus & Corsico (2004) about the presence of 
a diffusively-evolving double-layered chemical structure in 
pulsating DB white dwarfs. Our results show that clement 
diffusion proceeds very efficiently in these stars, causing 
the thickness of the pure helium envelope to increase by 
almost an order of magnitude by the time the DQ domain 
is reached, as compared with the situation at the begin- 
ning of the DB instability strip. 
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Table 1. Selected stages for our post-born again 0.5885 M0 sequence. 
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logTeff [K] 


log (4) 


logr [yr] 


log R, [cm] 


log g [cm s J 


log (xf) 




4.0466 


3.9940 


0.3629 


12.2700 


1.3524 


3.9908 


1.2555 


4.1044 


3.9928 


0.7537 


12.1539 


1.5848 


3.9884 


1.2555 


4.1493 


3.9916 


0.9283 


12.0634 


1.7657 


3.9865 


1.2555 


4.2061 


3.9898 


1.0999 


11.9489 


1.9948 


3.9836 


1.2554 


4.2586 


3.9878 


1.2323 


11.8429 


2.2067 


3.9804 


1.2554 


4.3060 


3.9856 


1.3402 


11.7470 


2.3985 


3.9770 


1.2554 


4.3586 


3.9826 


1.4526 


11.6402 


2.6120 


3.9724 


1.2554 


4.4053 


3.9795 


1.5480 


11.5453 


2.8019 


3.9675 


1.2553 


4.4552 


3.9756 


1.6473 


11.4435 


3.0056 


3.9611 


1.2553 


4.5000 


3.9714 


1.7348 


11.3519 


3.1888 


3.9542 


1.2549 


4.5537 


3.9655 


1.8392 


11.2416 


3.4092 


3.9442 


1.2548 


4.6017 


3.9593 


1.9345 


11.1425 


3.6075 


3.9328 


1.2548 


4.6500 


3.9519 


2.0355 


11.0421 


3.8083 


3.9200 


1.2547 


4.7014 


3.9424 


2.1563 


10.9346 


4.0233 


3.9003 


1.2548 


4.7515 


3.9315 


2.3064 


10.8290 


4.2345 


3.8734 


1.2549 


4.8010 


3.9221 


2.5381 


10.7251 


4.4422 


3.8384 


1.2559 


4.8508 


3.9156 


2.7410 


10.6225 


4.6476 


3.8360 


1.2585 


4.9003 


3.9007 


2.8811 


10.5159 


4.8607 


3.8366 


1.2603 


4.9500 


3.8804 


2.9975 


10.4064 


5.0797 


3.8240 


1.2621 


5.0002 


3.8549 


3.1382 


10.2933 


5.3060 


3.7937 


1.2649 


5.0503 


3.8224 


3.2890 


10.1768 


5.5390 


3.7610 


1.2693 


5.1005 


3.7774 


3.4503 


10.0538 


5.7848 


3.7121 


1.2758 


5.1502 


3.7134 


3.6242 


9.9224 


6.0478 


3.6410 


1.2868 


5.2001 


3.6081 


3.8226 


9.7700 


6.3525 


3.5196 


1.3038 


5.2466 


3.2645 


4.1517 


9.5053 


6.8819 


3.0875 


1.3515 


5.2001 


2.6917 


4.4184 


9.3118 


7.2688 


2.1510 


1.4052 


5.1507 


2.3614 


4.5701 


9.2454 


7.4018 


1.5060 


1.4294 


5.1001 


2.0635 


4.7480 


9.1978 


7.4970 


0.9067 


1.4455 


5.0525 


1.8007 


4.9492 


9.1615 


7.5694 


0.3861 


1.4480 


5.0018 


1.5620 


5.1430 


9.1435 


7.6055 


-0.2817 


1.3667 


4.9512 


1.2973 


5.4055 


9.1125 


7.6675 


-0.7823 


1.3114 


4.9011 


1.0430 


5.6464 


9.0855 


7.7215 


-1.2399 


1.1936 


4.8517 


0.7988 


5.8578 


9.0621 


7.7684 


-1.6699 


1.0408 


4.8003 


0.5493 


6.0509 


9.0403 


7.8120 


-2.2074 


0.8467 


4.7502 


0.3148 


6.2117 


9.0232 


7.8461 


-2.9096 


0.6406 


4.7011 


0.0919 


6.3536 


9.0100 


7.8726 


-3.7778 


0.4345 


4.6511 


-0.1312 


6.4955 


8.9984 


7.8957 


-4.8588 


0.2115 


4.6015 


-0.3477 


6.6278 


8.9893 


7.9139 


-6.0306 


-0.0137 


4.5540 


-0.5528 


6.7585 


8.9817 


7.9291 


-7.3201 


-0.2475 


4.5045 


-0.7646 


6.9069 


8.9748 


7.9429 


-8.9760 


-0.5207 


4.4558 


-0.9713 


7.0723 


8.9689 


7.9547 


-99.0000 


-0.8388 


4.4035 


-1.1917 


7.2858 


8.9633 


7.9660 


-99.0000 


-1.2737 


4.3540 


-1.3994 


7.5261 


8.9584 


7.9757 


-99.0000 


-1.8102 


4.3030 


-1.6132 


7.7839 


8.9536 


7.9853 


-99.0000 


-2.4829 


4.2505 


-1.8326 


8.0244 


8.9489 


7.9948 


-99.0000 


-3.2861 


4.2024 


-2.0322 


8.2128 


8.9454 


8.0018 


-99.0000 


-4.1214 


4.1527 


-2.2370 


8.3788 


8.9423 


8.0080 


-99.0000 


-4.9905 


4.1027 


-2.4419 


8.5257 


8.9398 


8.0128 


-99.0000 


-5.8260 


4.0500 


-2.6564 


8.6674 


8.9380 


8.0166 


-99.0000 


-7.9835 


4.0232 


-2.7651 


8.7349 


8.9373 


8.0180 


-99.0000 


-99.0000 



A novel aspect of our work is the fact that the calcu- 
lation of the evolutionary stages prior to the white-dwarf 
formation allows us to make sound predictions of the sur- 
face abundances expected in the DQ stars. In this connec- 
tion, we find that 12 C is by far the most abundant ele- 
ment that is convectively dredged-up to the outer layers, 



but with abundances far exceeding the low carbon abun- 
dances detected in many DQs (MacDonald et al. 1998) 
by a big margin. If we push our results to their limits 
(see below), then the plausibility of an evolutionary link 
between PG1159 and DQ stars with low detected carbon 
abundance appears to be unclear if convective dredge-up 
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Fig. 10. Abundances by mass of 4 He, 12 C and 16 as a 
function of the outer mass fraction for the 0.5885 Mq 
white-dwarf remnant having a final helium content of 
5.2 x 10 -3 Mq. Initial profiles at the start of the cooling 
branch (log T e g = 5.0) are plotted with thin lines. Later 
stages, denoted with normal and heavy lines, correspond 
to log T e ff = 4.41 and 4.28 respectively. 

is responsible for the observed carbon in DQs. Last but 
not least, our study also predicts the presence of trace 
amounts of 16 O, 13 C and 14 N in the atmospheres of these 
white dwarfs. 

However, it is conceivable that the occurrence of mass- 
loss episodes before and during the hot PG1159 stage 
could alter some aspects of the above mentioned results 
considerably. In this regard, recent observational evidence 
hints at the possibility that post-AGB mass-loss episodes 
could markedly reduce the mass of the helium-rich en- 
velope. Specifically, the immediate PG1159 predecessors, 
the [WC] stars, which roughly range between T g— 30000 
and 140000 K in the high-luminosity region of the early 
post-AGB evolution are known to have mass-loss rates 
of about 10 -55 to 10~ 6 - 5 Mq/vt that generally decrease 
with increasing T c ff. In addition, numerous hot, low- 
gravity PG1159 stars, have been reported (Koesterke & 
Werner 1998) to exhibit mass-loss rates ranging from 10 -7 
to 10 -8 Mq /yr. Finally, tentative evidence for the persis- 
tence of mass-loss rates of the same order along the hot 
end of the white dwarf cooling sequence and down to the 
domain of hot helium-rich white dwarfs has also been pre- 
sented (Werner 2001). Notably, the existence of PG1159 
stars with a helium content as low as 1 x 10~ 3 Mq has been 
suggested by asteroseismology in at least one of these stars 
with a stellar mass of 0.6 Mq (Kawaler & Bradley 1994), 
thus implying the occurrence of mass-loss during the evo- 
lution towards the PG1159 phase. These mass-loss rates 




Fig. 11. Logarithm of the number density of surface 12 C, 
16 0, 13 C and 14 N relative to that of 4 He as a function of 
the effective temperature. 



are indeed large enough to leave in principle signatures in 
the further evolution of these stars. 

To assess the possible implications of such mass-loss 
events for the white dwarf evolution, we have extended 
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Fig. 12. Same as Fig. 1101 but for a final helium content 
of 1.4 x 10~ 3 Mq. The initial profiles at the start of the 
cooling track (log T c g — 4.89) are plotted with thin lines. 
Later stages, denoted with normal and heavy lines, corre- 
spond to log T c ff = 4.407 and 4.29 respectively. 
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Fig. 13. Same as Fig. El but for a final helium content 
of 9.4 x 10~ 4 Mq. The initial profiles at the start of the 
cooling branch (log T c g = 4.93) are plotted with thin lines. 
Later stages, denoted with normal and heavy lines, corre- 
spond to log T c ff = 4.42 and 4.31 respectively. 

our calculations by considering two further evolutionary 
sequences in which mass loss is addressed during the 
hot post-AGB evolutionary stages. Specifically, we in- 
voke extreme constant mass-loss rates of 5 x 1CP 8 and 
1 x 10~ 7 A/o/yr along the hot post-AGB track from the 
low-gravity domain at T e ff« 115000 K sustained all the 
way down to the hot white dwarf cooling sequence at 
about T e ff=80000 K. As a result, we find that the he- 
lium content that eventually survives in the star amounts 
to 1.4 x 10~ 3 and 9.4 x 10~ 4 M Q , respectively. The total 
mass lost by the star amounts to about 0.013 Mq in both 
cases. Most of the mass is lost after the star has reached 
the point of maximum T ff. We note that mass loss un- 
covers deep regions in the star where the chemical com- 
position varies with depth. Thus, the outer layer chemical 
stratification after the end of mass loss at the start of 
the white dwarf cooling track looks somewhat different 
from the situation in which mass loss is not considered. In 
particular, the external chemical interfaces are markedly 
smoother. The implications for the chemical profiles ex- 
pected during the DB instability strip are clearly visual- 
ized with the help of Figs. El an d El f° r the mass-loss 
rates of 5 x 10~ 8 and 10 -7 Mq /yr, respectively. Note in 
particular that, in the case of the lowest helium content, 
a single and not a double-layered profile is more appro- 
priate to describe the outer layer chemical structure of 
pulsating DB white dwarfs. This is in line with the con- 
clusion arrived at in Althaus & Corsico (2004) that if post 
born-again DB white dwarf progenitors are formed with 
a helium content smaller than 10 -3 Mq a double-layered 
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T ff [10 3 K] 

Fig. 14. Number density of surface 12 C, relative to that of 
4 He as a function of the effective temperature. Solid and 
dotted lines correspond to the white dwarf remnant with 
helium content of 5.2 x 10~ 3 and 1.4 x 10~ 3 Mq, respec- 
tively. Filled circles refer to observed carbon abundance 
as taken from MacDonald et al. (1998). 

structure is not expected by the time the star reaches the 
red edge of the DB instability strip. Thus, the calculations 
presented here place that conclusion on a more solid basis. 
It is also clear that the initial chemical profile is different 
according to whether mass loss actually occurs or not, an 
aspect which is expected to bear some signature on the 
pulsational properties of variables DBs. 

We want to mention that, even in the case of ex- 
treme mass-loss rates of 10 _7 AfQ/yr , the helium enve- 
lope is not completely removed. However, the rate is high 
enough to erode any vestige of 1 N before the remnant 
reaches T cS w 88000 K. The detection of abundant 14 N in 
some coolest PG1159 stars (for instance in PG0122+200 
at Toff ~ 80000 K; see Dreizler & Heber 1998) makes the 
persistence of such extreme mass loss rather unlikely. 

Additionally, we investigate the consequences of mass- 
loss episodes during the hot post-AGB phase for the sur- 
face abundances expected in DQ stars. The predicted sur- 
face 12 C abundance together with the observed carbon 
abundance in DQ atmospheres are shown in Fig. 1141 Our 
results indicate that if DQ white dwarfs are the descen- 
dants of the post-born again PG1159 stars, then their 
surface carbon abundance is not expected to exhibit a 
marked dependence on the helium content with which the 
white dwarf is formed. In particular, note that the car- 
bon abundance far exceeds the low carbon abundances de- 
tected in numerous DQs. This prompts us to suggest that 
the DQs with low carbon abundance cannot be linked to 
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Fig. 15. Abundance by mass of 1 H, 4 He, 12 C and 16 as 
a function of the outer mass fraction for the 0.5885 Mq 
white dwarf remnant at T c g= 27000 K. Here, post-born 
again mass-loss episodes have been neglected. Thus, traces 
of hydrogen remain in the star at the beginning of the 
white dwarf evolution. Note that element diffusion leads 
to the formation of a DA white dwarf with a pure hydro- 
gen envelope mass of about 5 x 10~ 9 A/q. In this partic- 
ular simulation, we have neglected the diffusion of minor 
species such as 13 C and 14 N. 



the PG1159 stars if canonical convective dredge-up is the 
source of carbon for such DQs. Instead, they appear more 
likely to be good candidates for an evolutionary connec- 
tion that link them with stars that have somehow avoided 
the AGB thermally pulsing phase such as the hydrogen- 
poor AGB manque or the RCrB stars, with normal carbon 
abundances. 

Finally, we comment on the fact that high mass-loss 
rates like those suggested by observations have conse- 
quences for the evolutionary time-scales of the star. We 
find that for M = 5 x 10~ 8 and lO~ 7 M /yr, the star 
takes about 2.25 x 10 5 and 1.5 x 10 5 yr, respectively, to 
evolve from T cff = 117000 K down to T eS = 81000 K. This 
is markedly shorter than the time (4 x 10 5 yr) needed 
for the star to evolve across the same T e g interval in 
the absence of mass loss. This reduction in the evolu- 
tionary time-scales induced by mass-loss is expected to 
yield larger rates of period change in pulsating GW Vir 
stars, thus helping to partially overcome the discrepancy 
between theory and the observed value in the pulsating 
star PG1159-035 (Costa et al. 1999). 

In closing this section, we comment on the role played 
by the small amount of hydrogen that survives the proton 
burning during the last helium thermal pulse. As previ- 
ously mentioned, the hydrogen content eventually remain- 



ing in the star amounts to 1.3 x 10 -8 Mq, with traces of 
this element reaching layers as deep as 3 x 10 -4 Mq be- 
low the stellar surface. After the born-again phase, the 
total time spent by the remnant in the red-giant domain 
amounts to about 300-400 yr. Consequently, a mass-loss 
rate larger than 10~ 6 Mg/yr would be enough to remove, 
during this time interval, the last vestiges of hydrogen-rich 
material left in the star. It should be noted as well that 
the development of a surface lacking any hydrogen could 
also occur during the further PG1159 evolution as a re- 
sult of hydrogen burning and the persistence of a constant 
wind of about 10 -9 Mq /yr down to the hot white dwarf 
cooling branch. In the results presented thus far, we have 
assumed that it is indeed the actual course of events. From 
the opposite point of view, however, that is in the case of 
much weaker or less persistent mass-loss events, we would 
expect the formation of a DA white dwarf with a very thin 
hydrogen envelope as a result of element diffusion. This is 
indeed borne out by Fig. 1151 which illustrates the chem- 
ical abundance distribution by the time the white dwarf 
reaches T e g= 27000 K. The inner hydrogen has been dif- 
fused outwards and has formed a pure hydrogen envelope 
of mass 5 x 10~ 9 Mq, if the extreme situation of no mass 
loss is assumed, turning the white dwarf into one of the 
DA type. Thus, under these circunstances, it is conceiv- 
able that the born-again episode could also give rise to 
DA white dwarfs with very thin hydrogen envelopes. 

5. Conclusions 

In this paper we have studied some relevant aspects for 
the evolution of hydrogen-deficient white dwarfs (hereafter 
referred to as DB white dwarfs) by means of new evolu- 
tionary models based on a complete and self-consistent 
treatment of the evolutionary stages prior to the white 
dwarf formation. Specifically, we focused on DB white 
dwarfs, the progenitors of which have experienced a born- 
again episode, that is a very late helium thermal pulse on 
the early white-dwarf cooling branch after hydrogen burn- 
ing has almost ceased. The inclusion of a time-dependent 
scheme for the simultaneous treatment of nuclear evolu- 
tion and mixing processes due to convection, salt fingers 
and diffusive overshoot has allowed us to perform a de- 
tailed study of the abundance changes throughout all of 
the evolutionary phases, particularly during the thermally 
pulsing AGB state and the extremely short-lived phase of 
the born-again episode, for which the assumption of in- 
stantaneous mixing becomes completely inadequate. Our 
calculations made use the double-diffusive mixing length 
theory of convection for fluids with composition gradients 
(Grossman & Taam 1996). The study can thus be consid- 
ered as an assessment of its performance in low-mass stars. 
In particular, we have concentrated on the evolution of 
an initially 2.7 Mq star from the zero-age main sequence 
through the thermally pulsing and mass- loss phases to the 
white dwarf stage. 

As for the core helium burning and thermally pulsing 
AGB phases, our main results are: 
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— The inner carbon-oxygen profile exhibits a sharp vari- 
ation around M r w 0.33 Mq induced by mechanical 
overshoot. 

— Our models experience the third dredge-up during the 
AGB. Also, overshooting from below the short-lived 
helium-flash convection region at pulse peak causes 
the intershell region below the helium-rich buffer to 
be enriched with abundant oxygen, in agreement with 
Herwig (2000). 

— We find that pockets of 13 C and 14 N are formed at the 
base of the helium buffer after the end of the dredge- 
up phase. During the interpulse period, the 13 C-pockct 
is radiatively burnt, whilst the 14 N-pocket is engulfed 
by the helium-flash convection zone during the next 
thermal pulse. 

For the born-again evolution, the following results are 
worth mentioning: 

— The convection theory of Grossman & Taam (1996) 
appears to provide markedly shorter born-again time- 
scales than previous computations; this without invok- 
ing any reduction in the convective efficiency. However, 
they are still larger than the evolutionary time-scale of 
the born-again Sakurai object. This is a preliminary re- 
sult, however, which should be explored in more detail 
than attempted here. 

— After the born-again phase, we find the occurrence of 
a double-loop path in the HR diagram; that is, the 
star reaches giant dimensions for the second time af- 
ter the onset of the last helium thermal pulse and 
before finally returning to the white dwarf cooling 
track. The amount of hydrogen remaining in the star is 
1.3 x 10~ 8 Mq. Most of the hydrogen envelope burning 
occurs in about 1 month. 

— After the born-again episode, 4 Hc, 12 C and 16 are by 
far the dominant surface species with mass abundances 
of ( 4 He, 12 C, 16 0) = (0.306,0.376,0.228). Amongst the 
main remaining constituents 13 C, 14 N and 22 Ne are 
found with mass fractions of 4, 1.2 and 2.1 %, re- 
spectively. This is in agreement with surface abun- 
dance patterns observed in most hydrogen-deficient 
post-AGB stars such as PG1159 stars. The high 14 N 
surface abundance we find is also in line with that de- 
tected by Dreizler & Heber (1998) in some PG1159 
stars. 

— During the early PG1159 stage, about 3 x 10 -4 Mq of 
13 C is processed via the 13 C(a,n) 16 reaction. Also, 
the helium content in the star is reduced from 8.5 x 
10 -3 to 5.2 x 10~ 3 Mq as a result of helium burning. 

Additionally, the inclusion of the various processes re- 
sponsible for the element diffusion expected during the 
white dwarf regime has enabled us to extend the scope of 
our calculations to the domain of the helium-rich, carbon 
contaminated DQ white dwarfs, and more importantly, 
to assess the plausibility of the evolutionary connection 
PG1159-DB-DQ. In this regard, our main conclusions are: 



— Diffusion processes lead to the formation of a double- 
layered chemical structure during the DB evolution. 
In particular, by the time the domain of the pulsat- 
ing DBs is reached, above the helium-, carbon- and 
oxygen- rich intershell, there exists a pure helium layer 
of 2.3 x 10~ 6 M Q , which increases to 1.9 x 10~ 5 Mq 
when the star reaches the DQ domain. 

— As a result of convective dredge-up at low T e g, we find 
superficial 12 C with abundances far exceeding the low 
carbon abundances observed in many DQs. Our re- 
sults suggest that the DQs with observed low carbon 
abundances cannot be linked to the PG1159 stars if 
canonical convective dredge-up is the source of carbon 
for such DQs. 

In closing, we believe that our study would be rele- 
vant for pulsational applications since it provides a self- 
consistent description of the composition profile for both 
the core and the envelope. As recently shown by Metcalfe 
et al. (2003) this is an important aspect as far as astero- 
seismological inferences about pulsating DB white dwarfs 
are concerned. Additionally, our new models constitute a 
physically sound frame for answering some open questions 
about pulsating GW Vir stars, the hotter pulsating pre- 
decesors of variable DB stars (Gautschy 1997). We deem 
that our full evolutionary models would help to shed new 
light on the physical processes responsible for the driving 
of pulsations in these stars. Work in this direction is in 
progress. 

Detailed tabulations of our post-born again 
model are freely available at the following URL: 
http : //www . f caglp . unlp . edu . ar/evolgroup/ 
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